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' Abstract. We prove that the electron densities of electronic eigen 

Cl^i functions of atoms and molecules are smooth away from the nuclei. 
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1. Introduction and Statement of the Results. 



I We consider an A^-electron molecule with L fixed nuclei whose non- 

^ ■ relativistic Hamiltonian is given by 
O 

Q^. N f L 



Z,Z 



i^k 



o 

i^: + V - + V 

Dh' ^ \xi-xA ^ \Ri-Rk\' 

I ■ l<i<j<N ' * l<l<k<L ' ' 

^ '. where R = {Ri, R2, ■ ■ ■ , Rl) G K^^, Ri Rk for A; 7^ /, denote the 

positions of the L nuclei whose positive charges are given by Z = 
{Zi, Z2, ■ ■ ■ , Zl). The positions of the N electrons are denoted by 
(xi, X2, . . . , xat) G where Xj denotes the position of the j-th elec- 



j_j . tron in and A = XljLi is the 3iV-dimensional Laplacian. 



X 

H , 

. , The operator Hnl(II, Z) depends parametrically on R, Z, L, and 

the last term in corresponds to the internuclear repulsion which 

is just an additive term. It will play no role in this paper and we will 
hence neglect it, i.e. from now on we let 

H = Hjy iill, Z) — internuclear repulsion. 

The operator H is selfadjoint on L^(M^^) with operator domain 
V{H) = W'^'^iR^^) and quadratic form domain Q{H) = W^''^{M.^^), 
see e.g. Kato p. 
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We consider the eigenvalue problem 
(1.2) Hifj = (-A + V)ij = Eijj, 

where the potential V is given by 

^ ^ Zi 1 

^"^*^^ ^ I -7> . P, I ~^ I -7> . - I ' 

j=l i=l I 'I l<i<j<iV ' -J' 

Of course the eigenfunction ip and the associated eigenvalue E depend 
parametrically on R, Z, L, N. Since the potential V has singularities 
for Xj = Ri and for Xi = Xj one cannot expect classical solutions. 
In 1957 Kato showed that any local solution ip to (|1.2|) is locally 
Lipschitz, i.e. 

G C°,(M3^) and |V^| e L^,(M^^). 

Kato also characterized the behaviour of such a solution near points 
where two particles are close to each other (Cusp conditions). (Exten- 
sions of these results on the regularity of ip can be found in M. and 
T. HofFmann-Ostenhof and Stremnitzer ||^ and M. and T. Hoffmann- 
Ostenhof and 0stergaard S0rensen |Q). Of course, away from the sin- 
gularities of V, any local solution ip of ( p..2|) is smooth by elliptic reg- 
ularity. 

Equation (|1.2|) is a partial differential equation in 3A^ variables and 
hence only certain one-electron cases can be solved analytically. (One- 
electron atoms and diatomic one-electron molecules with equal nuclear 
charges) . 

Most of the information about bound states of atoms and molecules 
that Chemists and Physicists want to know is encoded in equation ( |1.2D . 
(Of course one can go beyond ( p..2[ ), for instance allow for nuclear 
motion, include relativistic corrections, etc.). Since the electrons are 
Fermions the physically relevant wavefunctions ip have to satisfy the 
Pauli principle. This amounts to requiring that ip transforms accord- 
ing to some specific irreducible representations of the symmetric group 
&N- Our results will include this. 

Already in the early times of quantum mechanics in the 1920 's and 
the 1930's various attempts were made to replace the full 3A^-dimensional 
equation ( |1.2| ) by simpler (usually non-linear) ones in 3 dimensions 
(Thomas- Fermi theory and Hartree-Fock theory). The relation of these 
approximations with the N-electron Schrodinger equation has been 



analysed in detail for big atoms, see e.g. Lieb and Simon |10|, Lieb 



and Lieb and Simon 11 



One important observation and motivation for the development of 
these and other approximation schemes was the insight that in order to 
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calculate the energy E or one- and two-electron operator expectation 
values there is no need for the full wave function ip{xi, X2, ■ ■ ■ , xn) but 
only for the one-electron density p{x),x G M^, the two-electron 
density p2{x,x'), {x,x') G and for the one-electron density 
matrix 7i(x,a;'), {x,x') G MP. These quantities are defined as fol- 
lows: Defining 

ijj{x) = ^{Xi, . . . ,Xj-i,X,Xj+i, . . .,Xn) 

and 

ipi^j^x, x) = ip{xi^ . . . , X, . . . , Xj-i, X , . . . , a; at), 

and assuming without loss of generality that is real valued (since the 
coefficients of H are real), the functions p, P2, and 71 are given by 

N 

I, J, 



:i.4) p{x) = j2 1 ^jd^. 



and 



l<i^j<N ■ 



N 



;i.6) = 



Here dxj, respectively dxi j, means integration over all variables except 



X 



31 



repectively Xi,Xj. Of course p{x) = 7i(x, x). 



More recently very successful approximation schemes have been de- 
velopped called Density Functional Theories (DFT). These schemes 
use some non-linear functionals in which only the one-electron density 
p occurs and lead to surprisingly good approximations to ground state 
energies and molecular geometries (see e.g. Eschrig [jl|). However, their 
relation to the full Schrodinger equation remains unclear. 

It is therefore surprising that the eigenfunction ip and in particular 
the electron density p itself, defined in ( p..4|) , have only rarely been the 
subject of mathematical analysis (see and references therein). 

Here we prove, based on recent work [0, natural smoothness results 
about the quantities defined in (pT^), (|1.5|), and (|1.6|). 

We do not assume anything about the eigenvalue E — in particular, 
it could be an embedded eigenvalue. The only assumption is that ip 
satisfies the following decay estimate: 

(1.7) |^(x)| < ce-^l''! for all x G M^^, 
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for some c, A > 0. 



Remark 1.1. Since ip is continuous, ( |1 . 7| ) is only an assumption on the 
behaviour near infinity. For references on the exponential decay of 
eigenfunctions, see e.g. Simon [|12 



Remark 1.2. We only assume the exponential decay for simplicity. Our 
proofs extend to the case where the eigenfunction ip decays faster than 
polynomially. 

The main result of this paper is the following theorem: 

Theorem 1.3. Let ip he an eigenfunction of H , satisfying ([l.Tj) . Let 
p,P2, and 7i be as defined in (|1.4|), (|1.5D, and ([L.6D, and define fur- 
thermore n{x) = p2{x,x). 

Let D = G M^} C M^ and define 

S = ... ,Rl}xR^)u (M^ X . . . , Rl}) C 

Then 

(1.8) p,neC°°{R'\{Ri,... ,Rl}) 
and 

(1.9) p2 e C~(R^\ (SUD)),7i e C~(M^\S). 

Furthermore, all the derivatives of p satisfy an exponential decay esti- 
mate near infinity: Let 7 G N'^ 6e a multi-index. Then for all e > 
there exists a constant c = 0(7, e) such that 

(1.10) 

\d2p{x)\ < ce-(^-^)l^'l for all \x\ > max{|i?i|, . . . , \Rl\} + 1. 

Similar results hold for n-electron densities and for n-electron density 
matrices. 

Remark 1.4. For simplicity we will only prove Theorem for atoms 
(i.e. for L = 1, i?i = 0) and only indicate the necessary modifications 
for the molecular case. Moreover, we only explicitly treat the density 
p, the proofs in the cases of n, p2, and 71 being essentially the same. 
Finally, it is enough to prove smoothness of each term J ijjj dxj in the 
sum i \1.4\) , and we therefore restrict ourselves to 

J -ipf dxi = J 'ip'^{x, X2, . . . ,xn) dx2 ■ ■ ■ dx^- 
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2. Additional regularity of eigenfunctions. 

We will need to know something about the regularity properties of ip 
in order to conclude that p is smooth. In order to study the regularity 
of ip it is convenient to work in spaces of functions which are Holder 
continuous. Let us recall the definition of Holder continuity: 

Definition 2.1. Let ilhe & domain in M", A; G N, and a G (0, 1]. We 
say that a function u belongs to C^'°'{Q) whenever u e C^(f2), and for 
all /3 e with \(3\ = k, and all open balls B{xo, r) with B{xo, r) C Q, 
we have 

sup J 1 _ |, <C(a;o,r). 

For any domain Q', with Q' C Q, we now define the following norms: 

Idf^uix) - d^u{y)\ 



= E ||9^u||L°°(n') + E sup 



m<k |^|^^x,yec',x^s, \x y\ 

We will need the following result on elliptic regularity in order to con- 
clude that the solutions of elliptic second order equations with bounded 
coefficients are C^'". The proposition is a reformulation of Corollary 
8.12 in Gilbarg and Trudinger [Q, adapted for our purposes: 

Proposition 2.2. Let Q be a bounded domain in and suppose u G 
W^''^{VL) is a weak solution of Au + ^"=1 bjDjU + Wu = g in Q, where 
bj,W,g e L°°{Q). Then u e C^'^^iQ) for all a G (0,1) and for any 
domain fl' , il' C we have 

\u\c^,<^(n') < C( sup \u\ + sup |5f|) 
n n 

for C = C(n, M, dist(^]^ dn)), with 

max {1, ||6j||L°°(c), || W^IU°°(c), ||fi'||L=°(n)} < M. 

j=l,... ,n 

Our regularity result on ip is the following: 

Lemma 2.3. Let ip be an N -electron atomic eigenfunction satisfying 
the decay estimate 7^ . Let P and Q be a partition of {1, . . . , A^}.- 

{1,... ,N} = PuQ, png = 0, p^0. 

Define xp as 

Xp = ^— y^^Xj e R^, 



P 



jeP 

and let T be any orthogonal transformation such that T(xi, ... , xat) 
{xp,x') with x' e R^^-^. 
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Let e > be given and define Up C as the open set: 



Up = < (xi, ... ,xn) e 



t,3N 



Xj\ > e for j G P, 



(2.1) \xj - Xk\ > e for j e P,k e qY 

Then 

d2^{tp o T*) e C^'\TUp) for all 7 G 
Furthermore, the following a priori estimate holds: 
o T*)ixp,x')\ + \d2(i; o T*)ixp,x')\ 



< Ce-^l(^^'^')l for all (xp, x') G TUf 



for some C = C{j) > 0. 



Remark 2.4. One could interpret f/p as a (very large) neighbourhood 
of a singularity x° of V where the electrons with coordinates Xj with 
j Q sit on the nucleus and the electrons with coordinates Xj with 
j ^ P sit on each other away from the nucleus, i.e. x*^ — (-^i; ' ' ' 1 -^a'')' 
with 



X 



x^j =0 for all j G Q, 
?/ V for all J G P. 



Notice that xp is (up to a scalar multiple) the centre of mass of the 
electrons j with j G P. Loosely speaking Lemma p.3| then says that in 
the neighbourhood Up, ip is smooth with respect to the centre of mass 
coordinate xp. 

One can also consider x*^ as a two-cluster singularity — one group of 
electrons on each other at the nucleus, another group of electrons on top 
of each other away from the nucleus. It is, of course, possible to have 
many clusters. Lemma \2,.'d\ is a special case of a more general theorem: 



If we have any number of clusters which are separated from each other 
and from the nucleus, then the eigenfunction ip can be differentiated 
any number of times with respect to the centre of mass of each cluster. 

We only need the two-cluster version (see Lemma p.3|) in the present 
paper. The more general result will be used in a detailed investigation 
of the regularity of the wavefunction near all kinds of singularities of 
the potential (see 0). 
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Remark 2.5. In the case of molecules the definition of Up is slightly 
different: 



MOL 



<>3N 



min — > e for j G P, 



\xj — Xfcl > e for j E P,k E q|. 
Apart from that, Lemma ^.3| remains unchanged. 



Before we prove Lemma let us fix some notation. We may as- 
sume without loss of generality that P = {!,... ,A^i}, with Ni < N. 
Then any orthogonal transformation T satisfying the assumptions in 



the statement of Lemma |2^ can be written as: 

1 



1 







\ 



V 



/ 



with the first row being understood as 3 x 3 matrices — first A^^i repeti- 
tions of -^7^ -^3 and then N — Ni repetitions of the 3x3 0-matrix. The 

remaining part of the matrix, T G M37v-3,3Ar(M) is such that the com- 
plete matrix T is orthogonal. We will denote the (3A^ — 3) x 3-columns 
of T by tj, i.e. 

f={h ■■■ t^), 
. Then we get: 

/ 



with tj G M3iv_3,3 



(2.2) 



\ 











^1 



%i+l 



i* 



/ 



Proof. For the proof of Lemma p.3| we first proceed as in M. and 
T. Hoffmann-Ostenhof and 0stergaard S0rensen We make the 
'Ansatz' 



(2.3) 



-,F-F, 
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with 



5Z 2'^^'"^ 5Z 4!^^' 

j=l l<j<k<N 



Xk\ 



and 



j=l l<j<k<N 



Observing (see with L = 1, /?i = 0), 

(2.4) AF = V, 

we get the following equation for tpi using ( |1.2| ) and (p73|): 

(2.5) At/^i + 2V(F - Fi) ■ VV^i + (|V(F - F^)\^ - AFi + F)^i = 0. 

Furthermore, Fi has the same behaviour as F at infinity, so we get the 
following estimates: 

FlWl^^^iN), \\F - Fl||ioo(]K3]V), ||V(F - Fl)||ioo(]R3]V) 

(2.6) < C(iV,Z, /?),/? G >0. 

We will first investigate the necessary regularity properties of the 
prefactor e^~^^ oT* in TUp with respect to xp. Then we will differen- 
tiate the equation (|2.5| ) with respect to and analyse the regularity 
of 92^(V'i oT*) in TUp. 

In the coordinates defined by T we have 

(2.7) X. = 4^ + t*x' for J < A^i, 

(2.8) Xj = t*x' forj>iVi. 
So 

— Xk = (t* — tl)x' for J, k < Ni or j, k > Ni. 

In particular these last expressions are independent of xp. So when we 
differentiate F o T* with respect to xp, the only non- vanishing terms 
come from derivatives of \xj\ o T* with j < A^i and \xj — Xk\ o T* 
with j < Ni, k > Ni. The definition of t/p implies that the function 
X 1-^ \xj\ is smooth with bounded derivatives on Up for j < A^^i, and 
5C I ^ I CC j CC I IS smooth with bounded derivatives on Up for j < Ni, 
k> Ni. Thus 

WVd^F - Fi) o T*\\l^^tu,) + \\dUF - Fi) o T*|U^(Tt;,) 

(2.9) < C(7) for all 7 G 
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Furthermore, we get using ( p.6|) that 

(2.10) 82^ {e^-^' o T*) G C°'\TUp) n L°^{TUp) for all 7 G 
Hence, due to ( |2.3| ) and (|2.9[), it remains to study the regularity of 

In the rest of the proof we will use the following notation: We will 
write F, Fi, ipi instead of F o T*, Fi o T*, tpi o T*. In particular we have 
the following relation 

Notice that since the Laplacian is invariant under orthogonal trans- 
formations we have the following equation for ip oT*: 

-A(V^ o T*) + {V o T*){ip o T*) = E{ip o T*). 

Now, once again using the invariance of the Laplacian and (p.4]) , 

AF = V oT\ 

so we get the following equation for ipi (compare with ( p.5|) ): 

L^Ji = 0, 

(2.11) L = A + 2V(F - Fi) • V + (|V(F - Fi)^ - AFi + E). 

The analysis of ipi will be based on the elliptic regularity result from 
Proposition p.2| . We shall proceed by induction and for this we will 
need the following open sets indexed by / G N: 

Ui = |(a;i, ... ,XN)e R^^ \xj\ > e(l - 2-('+i)) for j G P, 

\xj - Xk\ > e(l - 2-('+^)) for j EP,kEQy 

It is clear that for li < I2 we have 

Up C Ui, C Ui, C t/o. 
We will prove the following statement: 

Statement. For all 7 G N'^ we have 

(1) dl^^i G Ci'"(Tf/|^|) for all a G (0, 1). 

(2) dl^i,, G <','(T[/|,|). 

(3) There exists c = 0(7) > such that 

|V92,^i(xp,x')| + \dl^^i{xp,x')\ < ce-"l(^-^')l 
for all {xp,x') G T[/|^|. 
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Let US start by noticing that once this statement is estabhshed 
Lemma |2.3| is proved. 

The proof of the statement proceeds by induction with respect to 

In order to prove the statement for I7I = 0, let us look at the equation 
( |2.11| ). We get from Proposition |2.2| that 

since the coefficients of the equation are bounded on TUq, due 

to ( p.6|) . This proves (|ID for I7I = 0. 

We next prove that ipi G W^^^ (TUq). This is accomplished as follows: 
If we use that tpi G C^'°'{TUo) and that the coefficients in the equation 
( |2.11| ) (derivatives of F — Fi) are bounded (again using (^.6D), then it 
is easily seen from ( |2.11| ) that 

A^i e lUtUo). 

Therefrom, we get via standard elliptic regularity results (see for in- 
stance Folland @, Theorem 6.33]) that tpi G Wy^f{TUo). 

Next, we verify the exponential decay estimate (^) for I7I = 0. We 
know from the assumption (|1.7| ) that il^i decays exponentially, but we 
also need to prove it for Vipi- This is done exactly as in the induction 
step below, using the exponential decay estimate ([L^ ). In order not to 
repeat the argument, we refer the reader to the induction step below. 

Suppose now that we have proved (|l|)-(0) for all 7 with I7I < k. Take 
a 7 with length I7I = + 1. 

Differentiating the equation ( p.ll|) for ipi we get the following equa- 
tion for d^pipi (in the sense of distributions): 

L = A + 2V(F - Fi) ■ V + (I V(F - Fi)p - AFi + F), 

(2.12) + [d:^{\V{F - F,)\' - AFi + F)] d>ijr). 

In only partial derivatives of ipi of length < k occur. 

It is clear that in (|2.9|) we may replace TUp by TUq. Therefore we get, 
using the induction hypothesis for (|l|), that e L°° (TU\^\_i) . A priori 
dj ipi only satisfies the equation ( p.l2|) in the distributional sense. In 



order to apply Proposition we need that djpipi G W^^^{U\^\-i). 
However, this follows from the induction hypotheses for (0), and the 
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definition of Sobolev spaces. Now, due to ( |2.9| ), the coefficients of L 
are bounded, and therefore we get from Proposition ^]2| that d^pipi € 



Ci'"(T?7|^l_i) for all a e (0,1), proving (|). This together with (pj) 
implies that G Lf^^(TU\^\_i) . Therefore, we get (invoking once again 
that d^pipi e Ci'°(rf/|^|_i) for all a E (0, 1) and (p|)). 

As before, via a standard elliptic regularity theorem (e.g. Folland 0, 
Lemma 6.32]), we see that d2p'ijji G W^^^{TU\j\^i), and therefore (0) is 
proved. 

Finally, we prove the exponential decay of d].pipi- Let us write e; = 

e/2'+l Then we have for all x e f/z that 5(x,2e/) C f/z_i. 

Since the equation ( |2.12| ) is satisfied on T[/|^|_i it is in particular 
satisfied on TB{-x,2e\y\) for any x G f/i^y]. Applying Proposition |2^ 
with fi' = TB{:k, e|7|), = TB{^, 2e\^\), we get that 



l^xpV'ilci.-(n') < sup \d^piJi\ + sup 1/^ 

Notice, using ( |2.9|) , that the coefficients in the differential operator L 
are uniformly bounded on all of TUq. Furthermore, the derivatives of 
F and Fi in the expression for are also uniformly bounded on Uq — 
once again using (|2.9|). Therefore, using the induction hypothesis (^, 
we get for all {xp,x') G TU\^\: 



< C sup|a^^?/^i| +sup|/^ 

V o n 

< C sup e-^l(^^'^')l 

(ypy)eTB((xp,x'),2e|^|) 

In this last line the constant C does not depend on the position {xp, x') 
of the ball, since the coefficients of the equation ( |2.12| ) are uniformly 



bounded on all TUq. This proves (^ and thus finishes the induction. 

□ 

Remark 2.6. In the case of molecules we modify F and Fi as follows: 

^ ^ Zi 1 

1=1 j=l l<j<k<N 



12 S. FOURNAIS, M. AND T. HOFFMANN-OSTENHOF AND T. 0. S0RENSEN 

and 

Zi 



L N 

1=1 i=i 



Xj - Ri\ 



l<j<k<N 



4 v '^^ 



The rest of the proof is analogous. 



3. The proof of Theorem 1.3 



As noted in Remark |LJ it is enough to prove smoothness of each 
individual term in ( |1.4D . We therefore redefine p by 



xn) dx2 ■ ■ ■ dx 



TV- 



Lemma |2.3| will be essential in order to prove the smoothness of p. 
It suffices to prove that p e C°°(R3 \ S(0, i?)) for all R>0. Therefore, 
let us assume that |a;| > i? > 0. 

Remark 3.1. In the case of molecules we assume 

min \x — Ri\ > R> 0, 

1<1<L 

and prove that p E C^iM.^ \ (u[^-^B{Ri, R)^). 

Let Xij X2 be a partition of unity in M+i Xi + X2 = ^, Xi{^) = ^ on 
[0,R/{AN)], suppxi C [0,R/{2N)] and Xj e C-(M+) for j = 1,2. 
We combine the Xj's to make a partition of unity in M.^^ . Obviously: 

1= n (xi+x2){\xj - Xk\). 

l<j<k<N 

Multiplying out the above product, we get sums of products of xi's 
and X2's. We introduce the following index sets to control these sums: 
Define first 

M = {ij,k)e{l,...,Ny\j<k}, 

and let 

/CM, 
J = M\I. 

Now define, for each pair /, J as above. 
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Then we get 

1= n (Xi + X2)(|a;j -Xfcl) = ^ 0/(x), 

^<j<k<N ICM 

where the sum is over all subsets I C M. 
Therefore we obtain: 

pi^i) = / 4''^ {xi,X2, ... ,xn) dx2- ■■ dXN 

= iIj'^{xi,X2,..., XN)(f)l{xi, X2, . . . , Xn) dx2 ■ ■ ■ dXN 

ICM 

I CM 

We will prove smoothness of each of the individual p/'s, so let us pick 
and fix an arbitrary I C M. Our strategy is to associate to / a sub- 
set P of {1, ... , N} such that Lemma ^]3| is applicable, i.e. such that 
supp (pi C Up with Up defined in ( |2.1| ) with a suitable e > 0. 

Remark 3.2. To motivate the determination of P and hence the coor- 
dinate xp with respect to which we are allowed to differentiate, let us 
consider the following example: Let = 3 and / = {(1, 2), (2, 3)} and 
J = {(1,3)}. Then we have on the support of (f)j that l^i — X2\ < ^ 
(due to the suppxi) and since > R, \x2\ > Further, due to 
the suppx2, "we have \xi — x^l > ^. Suppose now, we would choose 
P = {1, 2} and Q = {3}, then, according to (pTTD, 

Up = {{xi, X2, X3)||xi| > e, \x2\ > e, \xi - x^] > e, \x2 - Xs\ > e}, 

for some e > 0. But then supp0/ <^ Up, since supp 0/ contains points 
with X2 = x^. On the other hand one easily checks that the choice 
P = {1,2, 3} is the right one. This example shows that we cannot just 
choose P to be {1} U {j < N\ (1, j) G /}. On the other hand P cannot 
be too big: For = 3, / = {(1,2)}, J = {(1, 3), (2, 3)} it is easily 
seen that with P = {1,2,3}, supp0/ contains points with X3 = and 
therefore supp 07- <^Up. 

Physically speaking, we divide the N electrons into 2 clusters. The 
electrons j with j E P define the 'maximal cluster' of electrons con- 
taining the electron 1. This will be done via an equivalence relation 
below. Note that the (three-dimensional) variable xp = ^J^^ Yljcp^j 

is (up to a scalar multiple) the centre (centre of mass) of the maximal 
cluster. 

Let ~ denote the equivalence relation on {1, . . . , N}"^ generated by 
/ and let P denote the equivalence class of 1. 
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Explicitly this means that j ~ A; if either j = k or there exists a 
sequence ji, ... ,ji with js G {1, . . . , N} for 1 < s < / and with js ^ jt 
for s 7^ t, such that 

(i) e / or e /, 

(ii) Us,js+i) e I or {js+i,js) G /, for 1 < s < / - 1, 

(iii) {ji,k) E I or {k,ji) E I. 
Clearly 

(3.1) l<N-2. 

Thus P = I}, Q = {1,... ,N}\ P. Notice that P ^ 0. In 

order to be able to apply Lemma |2.3| , we have to show that supp (pj C 
f/p (with a suitable choice of e in the definition of Up). 

Let j E P, then j ~ 1 and we can choose a sequence ji, . . . ,ji 
according to the above. Taking into account |xi| > R we have 

I CC j I I OC ]^ ^ j I * 

Further, with 

i-i 



\xi Xj\ < \xi Xjj^l + ^ ] ^is+il ~^ 



Xj\, 



s=l 



the length scale of the cut-off's and ( p.l|) we obtain that 

(3.2) 3 EP^ SUPP0/ C {(xi,... ,Xn) G M^^ | Ix^] > i?/4}. 

Remark 3.3. In the case of molecules we get 

j E P ^ supp0/ C {(xi, . . . ,Xn) G I min |xj - Ri\ > -R/4}. 

Furthermore, suppose j E P, k E Q, then it is clear that (j, k) E J or 
{k,j) E J (because if (j, k) E I, then 1 ~ j ~ and therefore k E P), 
and therefore: 
(3.3) 

j E P,k E Q ^ supp0/ C {(xi, . . . ,Xn) G M^^ | |Xj - Xk\ > R/{4:N)}. 

Using (p.2|) and (|3.3| ) we see that supp0/ C Up, with e = ^ in the 
definition (|2.1|) of Up. Hence we get from Lemma |2.3| that ip o T* is 
(infinitely often) differentiable with respect to the coordinate xp = 
— ^ J2jeP '^^ the support of (pi oT* . 

Denote gi = ip'^(j)i, and note that all partial derivatives of 0/ are 
bounded. We get from Lemma |2.3| that 

d2^{gioT*)EC''\R'''), 
(3.4) \d2^{gioT*){xp,x')\ < c,e~^l("^'"')l for all (xp,x') G R^^. 
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Now, we are ready to prove the smoothness of the electron density. 
We calculate (using the notation from (|2.2| ), and (|2.7| ) and ( p.8| )): 



gi{x, X2, ■ ■ ■ , Xjy) dx2 ■ ■ ■ dXN 

r-3 

gi{xi, . . . , xn)S{x — xi) dxidx2 ■ ■ ■ dxN 

I {gj o T*){xp, x')5{x y= — t^x') dxp dx' 

{gi o T*){^/n'i{x - tlx'),x') dx' 



(3.5) 

Using (|3.4| ) and Lebesgue integration theory, we obtain via the chain 
rule: 

dlPi{x) =82 {j {gio T*) (V^(x - tlx'),x') dx') 

= (v^)''' / {d2^{gioT*)){^/N;{x-tlx'),x')dx'. 



This proves that p is smooth away from the nucleus. 

The exponential decay of the derivatives of pi is a consequence of 
(|3.4| ). This can be seen by a similar calculation as in ( |3.5| ) but in 
reversed order: 



\d2Pi{x)\ < (v^)'"'y |(52,(^7,oT*))(V^(a:-t*x'),a:') 
< c [ e-^\(^^(--''^-')'-'^\ dx' 



dx' 



j (e-^l"l or*)(xp,x')5(x - ^ dxpdx' 



Let us write (x, X2, ■ ■ ■ , Xn) = {x, z). Then for all e G (0, 1): 

\{x,z)\ = (l-e)|(x,z)| +e|(x,z)| > (1 - e)\x\ + e\z\. 
Therefore 

'"e-"l(^'^)lrf^<e-^(i-^)l^-il / e-^'\'\dz 



ce 



-X{l-e)\x,\ 



This verifies inequality ( |1.10|) and finishes the proof of Theorem |1.3| . □ 
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